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Abstract 

We present a detailed analysis of Ramanujan's most accurate approximation to 
the perimeter of an ellipse. 
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1 Introduction 

Let a and b be the semi- major and semi-minor axes of an ellipse with perimeter p and whose 
eccentricity is k. The final sentence of Ramanujan's famous paper Modular Equations and 
Approximations to n, jS], says: 

" The following approximation for p [was] obtained empirically: 




i 



where e is about 



3ak 



20 



68719476736 

Ramanujan never explained his "empirical" method of obtaining this approximation, 
nor ever subsequently returned to this approximation, neither in his published work, nor in 
his Notebooks Indeed, although the Notebooks does contain the above approximation 
(see Entry 3 of Chapter XVIII) the statement there does not even mention his asymptotic 
error estimate stated above. 

Twenty years later Watson jH] claimed to have proven that Ramanujan's approxima- 
tion is in defect, but he never published his proof. 

In 1978, we established the following optimal version of Ramanujan's approximation: 

THEOREM 1. (Ramanujan's Approximation Theorem) Ramanujan 7 s approxima- 
tive perimeter 



p R := 7r <(a + b) + 



3(a-bf 



10(a + b) + Va 2 + 14afe + b 2 



underestimates the true perimeter, p, by 



n(a + b) ■ 0(A) • A 



10 



where 



A := 



a — b 
a + b 1 



1.2) 



;i.3) 



;i-4) 



and where the function 9(X) grows monotonically in ^ A ^ 1 while at the same time it 
satisfies the optimal inequalities 



3 .,.,14/22 

— < (9(A) ^ 7T 

2 17 v ; 11 I 7 



1.5) 



□ 

Please take note of the striking form of the sharp upper bound since it involves the 
/22 \ 

number I — — n J which measures the accuracy of ARCHIMEDES ' famous approximation, 
22 

— , to the transcendental numbers! 
7 

COROLLARY 1. The error in defect, e, as a function of X, grows monotonically for 
< A 1. 

□ 

COROLLARY 2. The error in defect, e, as a function of the eccentricity, e, is given by 

, , , 1.6) 



e(e) := a < 5(e) 



19 



,20 
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Moreover, e(e) grows monotonically with e, ^ e ^ 1, while 5(e) satisfies the optimal 
inequalities 



This Corollary 2 explains the significance of Ramanujan's own error estimate in (1.1). 
The latter is an asymptotic lower bound for e(e) but it is not the optimal one. That is given 



We sent an (updated) copy of our 1978 preprint to Professor Bruce Berndt in 1988 and 
he subsequently quoted its conclusions in his edition of Volume 3 of the Notebooks (see p. 
150 0). However the details of our proofs have never been published and so we have decided 
to present them in this paper. 

Berndt's discussion of Ramanujan's approximation includes Almkvist's very plau- 
sible suggestion that Ramanujan's "empirical process" was to develop a continued fraction 
expansion of Ivory's infinite series for the perimeter (P) as well as a proof, due indepen- 
dently to Almkvist and Askey, of our fundamental lemma (see §3). However, their proof 
is different from ours. 

The most recent work on the subject has been carried out by R. Barnard, K. Pearce, 
and K. Richards in |2] and was published in the year 2000. They also prove the major 
conclusion in our fundamental lemma, but their methods too are quite different from ours. 

3 Fundamental Lemma 

THEOREM 2. (Fundamental Lemma) Define the functions A(x) and B(x) and the 

coefficients A n and B n by: 




(1.7) 



□ 



in (1.7). 



2 Later History 



A(x) : 



1 + 



10 + V4 - 3x ' 



1 + A x x + A 2 x 2 H 



(3.1) 




x 



l + B x x + B 2 x 2 + ■■■ . 



(3.2) 



Then: 



A 1 = B l , A 2 = B 2 , 
A 5 < B 5 , A 6 < B 6 , 




(3.3) 
(3.4) 



where the strict inequalities in (3.4) are valid for all n ^ 5. 
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Proof. First we prove (3.3). We read this off directly from the numerical values of the 
expansion: 

1 



1 

16 
1 

64 
25 

4096' 



A 1 


= B 1 


A 2 


= B 2 


A 3 


= B 3 


A 4 


= B 4 



Now we prove (3.4) For A 5: B 5 , A e , and B 6 we verify (3.4) directly from their explicit 
numerical values. Namely, 



47- 

Ak = 4* B, 



49 

214 ' " u fjll' 

. 803 882 
A<5 = T^T' Bq = 



A 5 -B 5 = -±<0 
-79 

A 6 -B 6 = — <0. 



22i ' D 2 21 ' 
Therefore it is sufficient to prove 

A n < B n 

for all 

n > 7. 

Now the explicit formula for A n is 

A n = a n _i + a n _2 + a n -3 + ' ' ' + °i + a o 

where 



1 /2n-2 



a n-l 
On- 2 



2n - 3 16 n V n - 1 
1 1 /2ra - 2 



in— 1 



gn-2 | ~ 

2n - 5 16"- 1 V n - 1 J V 2 5 



ai 
a 



_J L( 2 W 2 M 

2 • 1 — 1 16 2 VI/ V 25 
16 V 2 5 



n-2 



Next we write 



and assert: 



/ -, . a n-2 . a n-3 . O-n-4 . 0\ Oq 

A„ = a n _i 1 H 1 1 1 1 1 



o n -i o n _i a n -i 



On-l 0>n-l 



(3.5) 
(3.6) 

(3.7) 



(3.8) 



(3.9) 
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CLAIM 1. 



The ratios a ' a k 1 decrease monotonically in absolute value as k increases from k = 1 

Q"n— k 

to k = n — 1. 

Proof. For k = 1, • • • ,n — 2, 

a n-k-l i | 



Qn—k 



2n-2k-3 



1 + 



An - Ak - 2 7 12 



^ - (which is the worst case and occurs when k = n — 2) 
< 1 



ai 



< 1. 



For k = n — I, 

This completes the proof. 
CLAIM 2. 



T/ie raizos fc 1 alternate in sign. 
Proof. This is a consequence of the definition of the a%. 



By CLAIM 1. and CLAIM 2. we can write (3.9) in the form 

A n = a„_i(l — something positive and smaller than 1) 
< a n -i. 

Therefore, to prove (3.8) for n ^ 7, it suffices to prove 
for all n ^ 7. 

By (3.8) and the definition of B n , this last afirmation is equivalent to proving 



□ 



□ 



(3.10) 



1 (2n-2 



gn-l 



< 



1 1 (In 



2n - 3 16 n V n - 1 / [ 2n - 1 4™ \ n 

which, after some algebra, reduces to proving the implication 

n 2n — 1 



n 7 >- 2 ^ n r 3 • 3" ' < I . 
2n 

n 



If we define for all integers n ^ 7 

n 2n — 1 



f[ii) . = 2 2n-3 . 3 n-i (3U) 



n 

then the affirmation (3.10) turns out to be equivalent to 

n^7 => f(n) <1 (3.12) 
This latter affirmation is a consequence of the following two conditions: 

1. /(7) < 1. 

2. /(7) > /(8) > /(9) > • • • > /(fc) > /(* + 1) > • • • 
Proof of 1. By direct numerical computation 

□ 

Proof of 2. We must show 

k>7^ f(k)> f(k+l). 

If we define 

**) == J^J-y (3.13) 

then we must show 

k^7=>g(k)>l. (3.14) 

Using the definition (3.11) of f(n) and the definition (3.14) of g(n), and reducing algebraically 
we find 



9(k) 



2k /2A; - l x 2 



6k -9 \ k+ 1 
and we must show that 

2 

6k -9 V k+ 1 
Define the rational function of the real variable x: 

2 



k^l^-^—C^-f) >1. (3.15) 
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Then the graph of y = g(x) has a vertical asymptote at x — | and 

lim g(x) = +00. (3-17) 



x 2 



Moreover, the derivative of ^(x) is given by: 

2(2x 2 - 7x + 1) 



9\x) 



x(x + l)(2x - l)(2x + 3) 



which implies that 



r<o if|<x<^ ; 

</(*)< = if* = 2*^, 
[>0 iix> 7 -±^. 

Therefore, for x > § g(x) decreases from "+00" at x = | (see (3.17)) to an absolute 
minimum value (in | ^ a; < 00) 

' 7 +^T) =1+ 37-y/iI 1 0363895208 ,.. 
4 y 399 + 69^41 

and then increases monotonically as x — > 00 to its asymptotic limit y = |. Therefore 

> 1.03638- • -for all x > -,x ^ * 

=^ g(x) > 1.03638 • • ■ for all integers n ^ 2. 

=>> /(n) > (1.03638 • • • )f(n + 1) for all integers n ^ 2. 

=>- /(n) > /(ra + 1) for all integers n ^ 2. 

=► /(7) > /(8) > /(9) > • • • 

which implies that the condition 2. holds. Moreover we conclude that 

f(n) < 1 for all integers n ^ 7 
=^(3.10)holds for all integers n ^ 7 
=^(3.5) holds for all integers n ^ 7 
=K3.4) holds for all integers n ^ 5 

and this completes the proof of the Fundamental Lemma. □ 
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4 Ivory's Identity 

In 1796, J. Ivory ^] published the following identity (in somewhat different notation): 

THEOREM 3. (Ivory's Identity) If ^ x ^ 1 then the following formula for ~B(x) is 
valid: 



TV 



1 + 2y/x cos(20) + X (1(f) = ^2 



n=0 



1 1 /2n 



2n - 1 4™ V n 



x n = B(x) 



(4.1) 



Proof. We sketch his elegant proof. 



7T 



1 + 2v/xcos(20) + x 64 



-1 /*7T 00 



m=0 

oo 



2m-14 m U y ' 



m ^2nimtf> 



1 + y/xe 2 ^^ 1 + v^e" 2 ^ 
1 1 (2n\ 



£ 

n=0 



2n - 1 4 n V n 



n 2nind 



iS{^Cr)^}S{^Cr)^}jT- 

f 1 1/2 

^ 1 2n - 1 4 n V ! 

n=0 ^ v 



□ 

We will need the following evaluation in our investigation of the accuracy of Ramanu- 
JAN's approximation. 



COROLLARY 1. 

Proof. By Ivory's identity, 

B(l) 



(4.2) 



7T JO 
1 



7T JO 
4 

7T 



1 + 2v / lcos(20) + 1 
^2 + 2 cos(2<£) d(j) 
v/4cos 2 (0) # 



□ 
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5 The Accuracy Lemma 



THEOREM 4. (Accuracy Lemma) For ^ x ^ \, the function 

3x 



A(x) := 1 + 



10 + V4 - 3x 



underestimates the function 



n=0 ^ \ / J 



X' 



by a discrepancy, A(x) which is never more than 



more than —=x h : 
2 17 



4 14 

7T ~ IT 



2 17 



x 5 < A(x) ^ 



4 14 

7T ~ IT 



X' 



(5.1) 



(5.2) 



x 5 and which is always 



(5.3) 



4 14 



7T 11 



Moreover, the constants and — — x 5 are the best possible. 



2 17 



Proof. By the definition of A(x) and B(x) given in Theorem 1., the discrepancy A(x) is 
given by the series 

A(x) := B(x) - A(x) 

= (B 5 - A 5 )x 5 + (B 6 - A 6 )x 6 + • • • 
:= 5 5 x 5 + <5 6 x 6 + • • • , 



where, again by Theorem 1., 



On the one hand 



5k > for k — 5, 6, • • • . 

A(x) =x 5 (5 5 + 5 6 x + ---) 

^x 5 (5 5 + 5 6 + 5 7 + ---) 
= x 5 A(l) 

= x 5 {B(l)-A(l)} 



x 



4 14 

7T ~ IT 



where we used Corollary 1 of Ivory's identity in the last equality. Therefore 







A(x) < (i - 
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/4 14\ 

This is half of the accuracy lemma. Moreover the constant is assumed for x — 1 

\tt 11) 

and thus cannot be replaced by anything smaller, i.e., it is the best possible constant. 
On the other hand, we can write 



A(x) = x 5 {5 5 + G(x)}, 



where 



G(x) := SqX + S 7 x 2 + 



G(x) > for all < x ^ 1 , 
G(x) ^ as x -> 0. 



This shows that 



and that 



A(x) > 5 5 x 5 = —x 5 



x->o x 5 2 17 



This proves both the other inequality in the theorem and the optimality of the constant 
3 

#5 = — , i.e., that it cannot be replaced by any larger constant. 



This completes the proof of the Accuracy Lemma. 



□ 



6 The Accuracy of Ramanujan's Approximation 

Now we can achieve the main goal of this paper, namely to prove Ramanujan's Approx- 
imation Theorem. 

First we express the perimeter of an ellipse and Ramanujan's approximative perimeter 
in terms of the functions A(x) and B(x). 

THEOREM 5. If p is the perimeter of an ellipse with semimajor axes a and b, and if p R 
is Ramanujan 's approximative perimeter, then: 



p = 7r(a + b) • B 



PR 



7i(a + b) ■ A 



a + b 



a + b 



(6.1) 
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Proof. We begin with Ivory's Identity (§4) and in it we substitute x :- 
the integral becomes 



a — b 
a + b 



Then 



i r 



TV 



\ 



1 + 24 



a — b 
a + b 



cos(20) + 



a — b 



4 H 



a + b J n(a + 6) J 



[a 2 sin 2 6 + b 2 cos 2 . 



and therefore 



B 



a — b 



a + b J J 7r(a + b) J 



(a 2 sin 2 d> + b 2 cos 2 > 



But, it is well known (BERNDT [3]) that the perimeter, p, of an ellipse with semiaxes 
a and b is given by 



p = 4 / (a sin 4> + b cos 

and thus 



" 2 • 2 , , .2 „ 2 







a — b 



p = n(a + b)-B{ ( — ) |>. (6.2) 



Moreover, some algebra shows us that 



a — o \ | \ a + o 



i,\ 2 
a — b 



10 + a/4-3 



a + 6 



1 f. 3(a-6) 2 
J (a + b) + 



a + 6 1 10(o + 6) + y/a? + Uab + b 2 

and we conclude that Ramanujan's approximative formula, pa is given by 

!>,< = x{a + b)A{ }■ (6.3) 



a + b 



The formula for p above was the object of Ivory's original paper 
Now we complete the proof of Theorem 1. 

Proof. Writing 

a — b 
11 



□ 



and using the notation of the statement of Theorem 1. we conclude that 



€ 



n(a + b) ■ 9(X) ■ A 



10 



where 




(6.4) 



□ 
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